TDA Progress Report 42-111 


November 15, 1992 






Mathematical Modeling of Bent-Axis Hydraulic 

Piston Motors 

R, D. Bartos 

Ground Antennas and Facilities Engineering Section 


Each of the DSN 70-m antennas uses 16 bent-axis hydraulic piston motors as 
part of the antenna drive system. On each of the two antenna axes, four motors are 
used to drive the antenna and four motors provide counter torque to remove the 
backlash in the antenna drive train. This article presents a mathematical model 
for bent-axis hydraulic piston motors. The model was developed to understand the t 
influence of the hydraulic motors on the performance of the DSN 70-m antennas’ 
servo control system. 


I. Introduction 

Hydraulic motors are often selected to drive mechani- 
cal devices because they are able to provide large output 
torque over a wide range of speeds while being relatively 
compact and lightweight. When a hydraulic motor is se- 
lected to drive a device under precision control conditions, 
it is necessary to understand how a particular motor will 
affect the overall system performance so that system com- 
ponents and control algorithms can be chosen to achieve 
optimum performance. The effect of a particular motor on 
system performance is best evaluated through computer 
simulations. This article develops a mathematical model, 
which can be implemented within any standard computer 
simulation package, for a bent-axis hydraulic piston motor. 
A glossary is provided in the Appendix. 


II. Theoretical Model Development 

The configuration of a typical bent-axis hydraulic pis- 
ton motor has been provided in Figs. 1 [1], 2 [1], 3, and 4. 


This section develops mathematical models which describe 
the output torque and fluid flow rate through a bent-axis 
hydraulic piston motor. 

A. Hydraulic Motor Torque 

Based upon the geometrical relationships indicated in 
Fig. 3 and Fig. 4, the instantaneous torque generated by 
the fluid pressure forces is found to be given by 

n 

T p = ^ ARPi sin y cos 9{ (1) 

t=i 

where 


27r(i - 1) 

6i = 0i + — 

( 2 ) 

n 


= P A if cos 6i > 0 

( 3 ) 
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Pi =: Pb if COS 9i < 0 


(4) 


nO 1 


if n = 4k - 2 (9) 


fc = 0,±1, ±2, ±3, . . . (5) 

The variables used in Eqs. (1) through (5) are defined as 

A = area of one piston 
n = the number of pistons 
P A = pressure at motor port A 
Pb — pressure at motor port B 
Pi = pressure acting on piston i 
R = pitch circle radius of the cylinder block bores 
7 = motor housing angle 

Oi = angle of the ith piston with respect to the valve 
plate, as shown in Fig. 4 

Note that Eq. (1) only expresses the torque generated by 
the pressure forces and does not model the net torque out- 
put of the motor since losses due to viscous forces and 
Coulomb friction are not considered. Substituting the mo- 
tor displacement per shaft revolution given by 

D = 2nAR sin j (6) 


Tp ** Cl + ° 2 sin ("Y") if n = Ak ( 10 ) 


A P = P a -Pb (11) 

*=1,2,3,... (12) 

The nondimensional constants C\ and C 2 are computed by 
means of least-squares estimation using Eq. (7) to obtain 
the actual motor torque. Parameter values are presented 
in Table 1 through Table 3, along with normalized stan- 
dard deviation of the error between the approximations 
given in Eqs. (8) through (12) and the actual values given 
by Eq. (7). The estimation of the parameters was per- 
formed using 500 data points evenly spaced over a single 
ripple. The constant C 2 physically represents the fraction 
of the mean torque equivalent to the peak to peak ampli- 
tude of the ripple because 


f 2w 

/ [Ci + C 2 |cos(n0i)|] d 6 1 = 1 

Jo 


C1 + C2 cos (y-) do 1 = 1 (14) 


into Eq. (1) yields 


Tp = Pi cos Oi 


J Ci + Ci sin (~^') d9i = 1 (15) 


mean torque 


DAP 


where fl. and Pi are as defined in Eqs. ( 2) through (5). 


the pressure. forcffU&t 


plate used to regulate the fluid pressure in the valve block 
bores. Figures 1 and 2 illustrate the interaction of the 
valve plate and cylinder block. The shear stress of the 
fluid in the clearance between the cylinder block and the 
valve plate is given by 


where 

h « r (23) 

and where the variables are defined as 


T = 


fjujr 


(18) 


where the variables are defined as 

h = clearance between the cylinder block and the 
valve plate 

r = radius from the center of the cylinder block 
H = absolute fluid viscosity 
w = angular velocity of the cylinder block 


h = clearance 

L = length of the cylinder block 
r = outside diameter of the cylinder block 

From Eqs. (21) and (22), it is evident that the torque 
due to viscous friction is proportional to both viscosity and 
angular velocity. Hence, the torque due to viscous friction 
can be expressed as 


Ty — — Ky DfiU) 


(24) 


The torque due to viscous forces is found by integrating 
the shear stress over the area of the cylinder block to valve 
plate interface shown in Fig. 5 by using the equation 

T-- I / ( rr)(rdipdr ) (19) 

J R\ J ip i 

Combining Eqs. (18) and (19) gives 



which after integration yields 

r= -£> 2 -*)(«;-«:) (2i) 

4 h 

To compute the total viscous torque of the valve plate to 
cylinder block interface, the results gathered from Eq. (21) 
for all surfaces where viscous forces exist must be added 
together. Viscous forces also exist in the clearance between 
the cylinder block and the motor housing. The resisting 
torque created by these forces can be found by integrat- 
ing the shear stress given by Eq. (18) over the lateral area 
of the cylinder block. By approximating the lateral sur- 
face of the cylinder block as a cylinder and integrating the 
shear stress over the lateral area of the cylinder, the torque 
created by viscous forces is found to be 

j< _ 27rpLr 3 u, (22) 

h 


where Ky is a constant of proportionality for a given mo- 
tor. The motor displacement, D , is included as part of 
the proportionality constant in Eq. (24) for mathematical 
convenience when one is experimentally estimating model 
parameters. 

An examination of Fig. 3 reveals that the five bearings 
of a hydraulic motor are loaded by the pressure forces act- 
ing upon each of the motor pistons. A forces and moment 
analysis on the shaft and piston subassembly shows that 
the bearings will be loaded proportionally to the hydraulic 
pressures in motor ports A and B, assuming that the pres- 
sure in the motor case is negligible. Since the motor bear- 
ings are not ideal and have Coulomb friction, there exists 
an average frictional torque acting on the motor shaft ac- 
cording to the equation 

T F1 = -D(P a + Pb){Kf + I<FS |sgn(w)|)sgn(w) (25) 
where 

{ -1 x < 0 

0 x = 0 > (26) 

1 x > 0 J 

and where the variables are defined as 

K f = a proportionality constant under non-stick-slip 
conditions 

Kfs = a proportionality constant under stick-slip con- 
ditions 
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Shaft torques resulting from Coulomb friction are pro- 
duced from the internal oil seals and from the piston mo- 
tion within the barrel of the hydraulic motor. This internal 
friction is not related to the fluid pressures existing within 
the motor and is given by the relationship 

Tf 2 - -T fc (1 + K s |sgn(w)|)sgn(w) (27) 

where T F c is the constant friction torque at high speed 
and K s is the proportionality constant under stick-slip 
conditions. 

The torque output available to the hydraulic motor load 
is the sum of the torque effects due to the piston forces, mo- 
tor inertia, viscous damping, and Coulomb friction. The 
net motor output torque to the load is therefore given by 
the equation 

T l =T p + Tj +T v + T F i + T F2 (28) 


B. Hydraulic Motor Flow 

Hydraulic motors are actuated by a control element 
such as a valve or pump which provides fluid flow to the 
motor. Hence, consideration of the fluid flow rate through 
the hydraulic motor is important from a system-level anal- 
ysis point of view. By computing the rate at which vol- 
ume is swept by the pistons using the geometrical relation- 
ships shown in Fig. 3 and Fig. 4, the theoretical flow rate 
through the motor, ignoring leakage effects, is found to be 


n 


Q A - ^2 I< A iAR sin 7 cos 

i = 1 

(29) 

n 

Q b — K gjARsm 7 cos OjW 

i = 1 

(30) 

K Ai = 1 if cos 0i > 0 

(31) 


The variables are defined as 

A = area of one piston 
K A i = a constant 
Kb% — a constant 
Q a = flow rate into motor port A 
Q B = flow rate into motor port B 
R = pitch circle radius of the cylinder block bores 
y = motor housing angle 
c v = angular velocity of the motor shaft 

Substituting the motor displacement given by Eq. (6) 
into Eqs. (29) and (30) yields 

Qa = ¥t-Y,Km coafc (36) 

2n “ 

»= i 

Qb = 77— X/ ^ Bi cos ( 37 ) 

ln .= i 

It is clear that Eqs. (36) and (37) are difficult to imple- 
ment within a computer simulation package because they 
involve very intensive computations. These equations also 
do not give any indication of the magnitude of the flow 
ripple relative to the mean value. By simulating the fluid 
flow through the motors using Eqs. (36) and (37), it is 
evident that the fluid flow rate through the motor can be 
approximated by the relations [2] 

Qa^^t- [Ci + C 2 |cos(n<?i)|] if n = 2k + 1 (38) 

2i r 

? DoJ f f3_, nJ ^ ndl \0_Jf - 2. 


£ 

II 

O 

if cos 9i < 0 

(32) 






Kbi — 1 

if cos 0i < 0 

(33) 

Duj 

17 

Cl + C 2 sin | 

7n0]\ 
V 2 ) 

if n = Ak 

(40) 


Kb, = 0 


if cos#* > 0 


(34) 


Qb = -Qa 


(41) 



The nondimensional constants C\ and C 2 are computed 
by means of least-squares estimation using Eqs. (36) and 
(37) to obtain the actual motor flow rates. It turns out 
that the nondimensional constants have the same values 
as those for the torque expressed in Eqs. (8) through (10). 
As a result of Eqs. (13), (14), and (15), the constant C 2 
physically represents the fraction of the mean flow equiv- 
alent to the peak to peak amplitude of the ripple where 


6 = slot width 
h = slot thickness 
/ = slot length 
A P = differential pressure 
H = absolute fluid viscosity 


mean flow rate = — — (43) 

2?r 

An evaluation of the flow rate through a hydraulic mo- 
tor must consider the internal leakage of the motor. The 
leakage past the hydraulic motor pistons from the pres- 
sure and return ports to the case drain can be modeled as 
laminar flow in an annulus between a circular shaft and a 
cylinder. The leakage flow rate is given by the equation [3] 


Applying Eq. (46) to a differential slot yields 

2*- r/» 3 dP_ 

® ~ 12 n dr 


(47) 


where the negative sign is used because flow occurs in the 
opposite direction of dP/dr. By rearranging Eq. (47) and 
integrating, the following expressions are obtained: 


Q = 


7T rc 3 

6 fiL 



(44) 



6/iQ 

7r/i 3 r 


dr 


(48) 


where 


P 


irh 3 


In r + C 


(49) 


c << r 

The variables are defined as 

c = radial clearance 
e = eccentricity of the piston 


(45) 


Applying boundary conditions and solving for the constant 
of integration, C, one obtains the relationship 


{Pl-P2)*h* 

^ in (f ) 


(50) 


L = length of the leakage path 

A P = the pressure difference between the motor cham- 
ber and the case 

r = radius of the cylinder bore 


Since Eq. (50) gives the radial leakage over an angle 
of 2 ir radians while the motor leakage only occurs over an 
angle <f>, Eq. (50) is adjusted by a factor of 4>/2w to express 
the leakage from the motor port to the case drain as 


During the operation of a hydraulic motor, fluid leaks 
between the valve plate and the cylinder block. The leak- 
age occurs from the case drain as shown by the curved 
arrows in Fig. 6. Also shown in Fig. 6 are the variables for 
one of the four leakage paths from the motor ports to the 
case drain. The leakage flow rate is derived as follows [4]: 
The leakage flow rate through a finite slot is given by 


^ (Pi - P 2 ) 

where the variables are defined as 
P x = pressure at radius 1 


Q = 


APbh 3 

12 pi 


(46) 


P 2 = pressure at radius 2 
Ri = radius 1 


where the variables are defined as 


R 2 = radius 2 


(51) 
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(56) 


<j) = angular length of kidney-shaped motor port in 
radians 


Q = 


bh 3 

12/i/ 


AP 


Using Eqs. (44) and (51), it is clear that the internal leak- 
age from the port to the case is proportional to the pres- 
sure difference between the port and the case and inversely 
proportional to the fluid viscosity. Because the length of 
the leakage path found in Eq. (44) changes with the an- 
gular shaft position, the leakage proportionality constant 
will be a function of the motor shaft angle. The internal 
leakage variation with shaft angle has been observed by 
Johnson [5]. Hence, the internal leakage from the motor 
ports to the case drain is expressed as 


Qla 


K(6) L a(Pa-Pc) 

P 


(52) 


Qlb 


K (®)lb (P& P ° ) 
P 


(53) 


where 

I<(0) LA = constant of proportionality for leakage from 
port A to case 

K (0)lb = constant of proportionality for leakage from 
port B to case 


Since there are two identical cross port leakage paths, the 
total cross port leakage is expressed as 


Q = 


bh 3 

6fil 


AP 


(57) 


where the variables are defined as 

6 — width of leakage path 
h — clearance 
/ = length of leakage path 
A P = pressure difference across the motor 

Equation (57) shows that the cross port leakage is propor- 
tional to the pressure difference between the motor ports 
and inversely proportional to the viscosity of the fluid. The 
leakage from port A to port B can be expressed as 


Qab — 


Kab(Pa-Pb) 

P 


(58) 


Pc = case pressure 

Because the valve plate is symmetrical across the line di- 
viding the A and B motor chambers, it is expected that 

K(0) LA = I<{e-v) LB (54) 

where <p is a phase angle given by 


where Kab is a constant of proportionality determined 
through experiments. 

The net flow passing into port A, port B, and the case 
of a hydraulic motor is determined by summing the effects 
of the rate at which volume is swept by the pistons and 
the effects of internal leakage. The net flow to port A and 
port B is 



and n is the number of pistons. 

Internal hydraulic motor leakage exists between the two 
motor ports. The leakage flow between motor ports occurs 
across the valve plate to cylinder block sealing surface as 
illustrated in Fig. 7. Assume that the leakage path be- 
tween the motor ports is modeled as a rectangular leakage 
path. The leakage flow rate for a single cross port leakage 
path is given by [6] 


Qan ~ Qa T Qla + Qab ( 59 ) 

Qbn — Qb T Qlb ~ Qab (50) 

respectively, and the net leakage to the case is 

Qcn = Qla + Qlb ( 51 ) 

Positive flow is defined as flow into the motor port or case. 
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III. Conclusion 

This article has presented a mathematical model for a 
bent-axis hydraulic piston motor; the motor can be used to 
perform dynamic system simulations which predict system 
performance under various conditions. The model param- 
eters must be experimentally determined for a particular 
motor before the model can be used effectively because mo- 
tor manufacturers do not possess this type of information 


or are unwilling to provide it. A discussion of the exper- 
imental test procedures required to determine the model 
parameters is beyond the scope of this article but can be 
obtained: see Bartos, 1 Merrit [3], and the Hydraulic Fluid 
Power documents [7,8]. 

1 R. D. Bartos, Statement oj Work 33S-SW-1015: Hydraulic Motor 
and Relief Valve Package Experimental Evaluation (internal doc- 
uxnent), Jet Propulsion Laboratory, Pasadena, California, January 
7 , 1992 . 
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Table 1. Nondimenslonal parameter values for n s 2k + 1. 


Number 

of 

pistons 

Ci 

C 2 

2^(7 /DAP 
or 

2n<j/Du> 

3 

9.1331e-01 

1.3617e-01 

2.1080e-03 

5 

9.6932e-01 

4.8183e-02 

8.0199e-04 

7 

9.8442e-01 

2.4468e— 02 

4.1507e-04 

9 

9.9059e-01 

1.4773e-02 

2.5254e-04 

11 

9.9371e~01 

9.8797e-03 

1.6955e-04 

13 

9.9550e-01 

7.0697e-03 

1.2159e-04 

15 

9.9662e-01 

5.3083e-03 

9.1425e-05 

17 

9.9737e—01 

4.1318e-03 

7.1228e-05 

19 

9.9789e-01 

3.3072e-03 

5.7050e-05 


C\ = a dimensionless constant. 

C 2 — a dimensionless constant. 

a = standard deviation of the approximation error. 

2k(t/DAP = standard deviation of the approximation error 
normalized with respect to the average torque. 

2 ir<r/Du! = standard deviation of the approximation error nor- 
malized with respect to the average flow rate. 


Table 2. Nondimenslonal parameter values for n = 4k - 2. 


Number 

of 

pistons 

Ci 

C 2 

2ira / DAP 
or 

2it(j /Du 

2 

2.9905e-16 

1.5708e+00 

6.5742e-16 

6 

9.1331e-01 

1.3617e-01 

2.1080e-03 

10 

9.6932e-01 

4.8183e-02 

8.0199e-04 

14 

9.8442e-01 

2.4468e-02 

4.1507e-04 

18 

9.9059e-01 

1.4773e-02 

2.5254^04 


Ci = a dimensionless constant. 

C 2 ~ a dimensionless constant. 

a = standard deviation of the approximation error. 

2'ko/DAP — standard deviation of the approximation error 
normalized with respect to the average torque. 


Table 3. Nondimenslonal parameter values for n = 4k. 


Number 

of 

pistons 

Ci 

c 2 

2tt 0 / DAP 
or 

2 tt<t/ Du> 

4 

7.9791e-01 

3.1740e-01 

4.2129e-03 

8 

9.5176e— 01 

7.5756e— 02 

1.2421e-03 

12 

9.7874e-01 

3.3390e-02 

5.6651e-04 

16 

9.8808e-01 

1.8727e-02 

3,2147e-04 

20 

9.9238e-01 

1.1970e-02 

2.0657e-4)4 


C 1 = a dimensionless constant. 

C 2 — a dimensionless const amt. 
a — standard deviation of the approximation error. 

/DAP = standard deviation of the approximation error 
normalized with respect to the average torque. 


SHAFT BEARINGS 





A PISTON AREA 

AXIAL REACTION FORCE OF BEARINGS 
F BRj RADIAL REACTION FORCE OF BEARING / 
Pi PRESSURE ACTING ON PISTON / 

R PITCH CIRCLE RADIUS 


J F TORQUE DUE TO OIL SEAL AND 
PISTON FRICTION 

Tv TORQUE DUE TO VISCOUS FORCES 
Y MOTOR HOUSING ANGLE 
0 PISTON ANGLE FROM REFERENCE 
OF PISTON / 


Fig. 3. Shaft and piston subassembly of a bent-axis hydraulic piston motor. 



Fig. 5. One region of fluid shear stress acting between the valve 
plate to motor cylinder block interface. 
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Fig. 6. Leakage llow paths across the valve plate sealin 
surfaces to the motor case. 



Fig. 7. Leakage flow path across the valve plate 
sealing surfaces between the two motor ports. 
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Appendix 


A = area of one piston 
b = slot width, width of leakage path 
c = radial clearance 
C — a constant of integration 
Ci = a dimensionless constant 
C 2 = a dimensionless constant 
£) = motor displacement per revolution 
e — eccentricity of the piston 
h = clearance, slot thickness 
z = piston number 
J = motor inertia 
= integer 
K a i — a constant 
I<Bi — a constant 

A'(0)yi£ = a constant of proportionality related to the 
leakage between motor ports 

K f — a proportionality constant under kinetic fric- 
tion conditions 

Kps “ a proportionality constant under static friction 
conditions 

K(Q)la - a constant of proportionality for leakage from 
port A to case 

K(6)lb = a constant of proportionality for leakage from 
port B to case 

Ks = a proportionality constant under static friction 
conditions 

K v — a constant of proportionality related to viscous 
friction 

/ = slot length, length of leakage path 
L — cylinder block length 
n = number of pistons 
P — pressure 

P A = pressure of motor port A 
P B = pressure of motor port B 
P c — case pressure 

Pi = pressureactingonpiston i,pressureatboundary i 


A P = differential pressure across the motor 
Q a = theoretical fluid flow rate into motor port A 
Qab = leakage flow rate from port A to port B 
Q an = net fluid flow rate into motor port A 
Q b = theoretical fluid flow rate into motor port B 
Q bn = net fluid flow rate into motor port B 
Q cn = net leakage flow to the case drain 
Qla = internal leakage from port A to case 
Q bb — internal leakage from port B to case 
r = radius 

R = pitch circle radius of the cylinder block 
Ri = radius at boundary i 
T — torque 

T F c = constant friction torque at high speed independent 
of port pressures 

Tn = torque due to Coulomb friction originating from 
motor bearings 

7>2 = torque due to Coulomb friction originating from 
oil seals and piston motion 

Tj — torque to accelerate the inertial load 

Tl = torque delivered to the load 

Tp = torque generated by fluid pressure 

Ty = torque required to overcome viscous friction 

x = a variable 

a = angular acceleration of the motor shaft 

7 = motor housing angle 

Oi = angle of piston i 

fi =: absolute fluid viscosity 

p = fluid density 

r = shear stress 

(p = phase angle 

<f> — angular length of kidney-shaped port 
^ — angle 

<jj = angular velocity of the motor shaft 
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